Energy and contact of the one-dimensional Fermi polaron at zero and finite 

temperature 
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We use the T-matrix approach for studying highly polarized homogeneous Fermi gases in one 
dimension with repulsive or attractive contact interactions. Using this approach, we compute ground 
state energies and values for the contact parameter that show excellent agreement with exact and 
other numerical methods at zero temperature, even in the strongly interacting regime. Furthermore, 
we derive an exact expression for the value of the contact parameter in one dimension at zero 
temperature. The model is then extended and used for studying the temperature dependence of 
ground state energies and the contact parameter. 
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Strongly interacting Fermi gases are of interest in many 
fields of physics, such as atomic and nuclear physics, as 
well as the study of neutron stars and quantum chromo- 
dynamics. Furthermore, a better understanding of Fermi 
gases is essential for improving our description of elusive 
condensed matter phenomena, such as high-temperature 
superconductivity and (anti)ferromagnetism. The Fermi 
polaron is part of a class of phenomena knovifn as impu- 
rity problems, which occur in a wide range of physical 
disciplines. 

Recently it was shown in a series of papers by Tan 
[l|,l3], building on important early contributions provided 
by Olshanii and Dunjko [3|] , that an interacting Fermi gas 
with short-range interactions exhibits a universal decay 
of the momentum distribution n^ ~ C/fc* as A; — > oo. 
The quantity C is called the contact and is intimately 
connected to the thermodynamics of the many-body sys- 
tem. In general it depends on physical properties such as 
temperature, pressure and scattering length. However, 
the universal decay holds for any value of these quan- 
tities as long as the interactions are short-range (for a 
review on Tan's contact parameter, see [J|). This behav- 
ior persists even in one dimension [5[. The contact has 
also been probed experimentally in Fermi gases [6|, lZ| and 
a Bose- Einstein condensate [8] . Since the contact of a ho- 
mogeneous Fermi gas can be measured [9|, this provides 
an excellent test of theories. In particular, our work is 
relevant for a ID quantum wire geometry jlOf . 

The ground state energy of a polaron, a particle with 
generalized spin down interacting with a bath of spin 
up-particles (in experimental practice, the two different 
"spin" states are generally different hyperfine states), 
is predicted surprisingly well using a simple variational 
Ansatz due to Chevy [11[. This work was later clarified 
and extended to other dimensions 12l4l4 1 and applied to 
a lattice geometry [l5|. Another siinple approach gives 
reasonable results in all dimensions [16|. In ID it is well 
known that the usual Fermi liquid picture breaks down 
due to the Peierls instability, and the quasiparticle pic- 
ture fails due to the importance of collective excitations 



17|. However, in the case of polarons the disturbance of 
the ideal Fermi distribution of the atoms of the major- 
ity component is negligible in the thermodynamic limit 
and we expect the diagrammatic approach to be applica- 
ble. Such an approximative approach may seem super- 
fluous given the exact results available from the Bethe 
Ansatz (BA) solution, but the T-matrix approach is sim- 
pler, more physically intuitive, easily extended to finite 
temperature and can be generalized to higher dimensions, 
where no BA solution is available. 

We apply the non-self-consistent (NSC) T-matrix ap- 
proach to the problem of a Fermi polaron in ID, at 
zero as well as finite temperatures and for repulsive as 
well as attractive interactions. In doing so, we compute 
the ground state energy of the polaron and the contact 
density. Few exact solutions are available for fermionic 
many-body systems; we derive an exact expression for the 
contact density based on the BA solution, valid at zero 
temperature. We compare our results to a more elabo- 
rate semi-self-consistent scheme based on the Brueckner- 
Goldstone theory ijjj, il9|] , which has been widely used in 
the context of nuclear matter [20| , and (perhaps counter- 
intuitively) find that the self-consistent approach is less 
accurate than the NSC approach. 

We consider a system of free fermions in ID with a 
"spin down" particle interacting through a (5-function po- 
tential with a sea of "spin up" particles. To wit, we use an 
interparticle potential (we set Ti — 1) V{x) = — —5[x), 
where m is the mass of a particle and a is the ID scat- 
tering length. The external potential is set to zero. The 
system is thus described by the Hamiltonian: 
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+ fdx f dx'^\{x)^l{x')V{x - x')i;^{x')i^^{x), (1) 



where a is the spin index, fj,^ is the chemical potential 
and ^cr (x) destroys (creates) a particle with spin a at 
position X. 

The off-shell many-body T-matrix describes all possi- 



ble (relevant) scattering events involving at most a sin- 
gle particle-hole-pair. It can be obtained numerically 
from the two-body T-matrix, which is easy to calcu- 
late for a (5-function potential. In momentum space, this 
potential is given by a constant and the two-body T- 
matrix ro(S') can be obtained either from the Lippman- 
Schwinger equation [2l| or from a scattering calculation 
using the scattering amplitude f{k) = j^^^ (k is mo- 
mentum) [24I, with identical results. The many-body T- 
matrix r{S), which in the case of a (5-potential depends 
only on the center-of-mass 1 -|- 1-momenta (momentum 
and frequency), is then used according to the Galitskii 
integral approach in the ladder approximation to obtain 
the self-energy: 

np) = -^ j -^Gn{K)ns), (2) 

where K — {fc,w} is the generalized 1 + 1-momentum, 
Go{K) is the noninteracting Green's function for the ma- 
jority component, and S = ^{P + K) are the center-of- 
mass variables. The ground state energy is obtained by 
starting with a guess for the polaron energy Ep and then 

iterating equation ([2]) using Ep = Cp + 'E{p,Ep) (where 

2 
Cp = ^ is the kinetic energy) to obtain a stationary so- 
lution. The Green's functions entering the T-matrix in 
the NSC theory are undressed, hence the iterative loop 
in eq. ([2]) does not alter the T-matrix. As we shall see, 
improving upon the NSC theory is not easy and it is 
a better approximation than what may be expected at 
first glance. To generalize this approach to finite tem- 
peratures, consider the noninteracting Green's function: 
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(3) 



where 9 is the step function as usual, rj ^ Q^ and /cp 
is the Fermi momentum of the majority component. In 
the thermodynamic limit, we can simply replace the step 
function by the Fermi-Dirac distribution at finite tem- 
perature. 

We compute the polaron ground state energy i?p=o hi 
ID by numerically evaluating E at zero momentum over 
a large range of interaction strengths centered around 
infinitely strong interactions, i.e. k^a = (see Figure 
[1]). Our results show good agreement with exact results 
obtained from the BA 23|, [24 1 as well as numerical re- 
sults in the infinitely repulsive limit fcpo — >■ 0~ obtained 
from a variational Ansatz [l3|. In fact, in this limit the 
NSC theory yields E^sc/E^ « 1.24 (for fcpa = -IQ-^), 
which is closer to the exact result of Ee^/Ep — 1 than 
the value reported in [13| of Evur/Ep « 1.39. So while it 
was argued that in 3D the first order variational Ansatz 
of Chevy and the T-matrix approach are the same at zero 
temperature [12|, this does not appear to be the case in 
ID. We argue that this feature is related to the fact that 
in ID there is no ultraviolent divergence or a renormal- 
ization requirement, so that scattering of holes within 



the Fermi sphere is still relevant in ID. In the attractive 
regime, our numerics converge for fcpa ^ 0.5. At smaller 
(positive) values of the interaction strength our approach 
breaks down, possibly due to numerical instabilities. 
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FIG. 1. Energy of the ground state of the ID Fermi polaron 
as a function of the interaction strength kpa. Solid line shows 
the exact (BA) result, crosses (NSC) and circles (Brueckner- 
Goldstone theory) are numerical results from iterating eq. Q 
at various values of fcpa. Error estimates for the Brueckner 
theory results are related to poor convergence. 



We have also calculated the ground state energies using 
the Brueckner- Goldstone theory 1^, [13], which adds a 
self-consistency requirement to the T-matrix. The result 
is significantly worse than the NSC theory (see Figurc[T|). 



This lends credence to the idea put forth in [25[, where 
the authors argued that the self-consistent treatment per- 
forms poorly in highly spin-imbalanced systems because 
of the near-perfect destructive interference between dia- 
grams that are not included in the NSC treatment. Ex- 
tending to a self-consistent treatment only includes some 
of these diagrams, and the destructive interference is lost. 

At finite temperature (see Figure[2]), we find that in the 
high-temperature limit, the mean-field limit E^f/ E-p — 
—A/nk-pa is recovered, as expected. Furthermore, as 
|fcpa| is reduced, it appears that higher temperatures are 
needed to reach the crossover from the BA-regime to the 
mean-field regime. 

The Tan's contact is a central quantity connecting var- 
ious thermodynamic quantities (such as energy, pressure 
and entropy) in atomic gases with short range interac- 
tions. It can be obtained from the asymptotic momen- 
tum distribution n^ as C = limfe_j.oo k^rik and is the same 
for both the minority and the majority component. Since 
the effect of the impurity on the momentum distribution 
of the majority component is vanishingly small, we can 
use the Fermi-Dirac distribution for the majority com- 
ponent, even though this does not have a 1/k'^ tail. The 
many-body T-matrix can be used for calculating the con- 
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FIG. 2. (Color online) Energy of the ground state of the Fermi 
polaron as a function of inverse temperature Tp/T (Tp is the 
Fermi temperature). Top (green) line shows the repulsive 
polaron at fcpa = —1, bottom (red) line shows the attractive 
polaron at fcpa = 1. Inset; normalized energy E' = ^~_^™ 
for various values of fcpa = —0.5, —1, —2, —10, 1, where "mf" 
is the mean-field prediction and "exact" is the ground state 
polaron energy at T = for the considered scattering length. 
The bottom (red) line is for fcpa — —0.5. 



tact from a perturbative approach as 26l |27|: 

^p-t-iij z^rU' — 2 — 



"■t,fc, 



(4) 



where £4,^0 is the polaron ground state energy and e-f^k — 
,2 
2^. Alternatively, we can also compute the contact both 

numerically and exactly by evaluating the derivative of 

the ground state energy with respect to the scattering 

length, i.e. using the ID Tan adiabatic theorem [5|: 



dE 
da 



C 
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(5) 



where to^ — m/2 is the reduced mass. Based on this 
theorem, we can use the exact BA ground state energies 
to calculate the exact dimensionless contact density Ccx = 
C/{Lkp), where L — ir/k-p^i is the length of the one- 
dimensional gas and k-p.i is the Fermi momentum of the 
minority component. The result is: 
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(6) 



In the infinitely repulsive limit, fcpa — > 0^, this ap- 
proaches the finite value of 4/37r^, which compares to the 
value 0.749... for the spin-balanced case [5|, and is exactly 
the same as the value obtained for a Tonks- Girardeau 
(TG) gas of infinitely repulsive one-dimensional bosons 
[28| . Thus, an infinitely strongly repulsive impurity has 
the same contact when interacting with a "fermionized" 
sea of bosons as with an ideal Fermi sea. This is a some- 
what surprising result considering that the momentum 



distributions of the TG Bose gas and the ideal Fermi gas 
are quite different, but one should observe that the im- 
purity and the ideal Fermi gas actually constitute a TG 
Fermi gas |29| . 
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FIG. 3. Inverse contact density 1/C as a function of interac- 
tion strength at zero temperature. Solid line: exact result ((6|. 
Dotted line: mean-field prediction. Crosses are numerical re- 
sults obtained using eq. (|4]), squares are obtained from a mod- 
ified version of eq. Q (see main text) and stars are obtained 
from an implementation of the ID Tan adiabatic theorem ((S)) 
using numerically obtained values for the polaron energy. In- 
set: exact and mean-field values of the contact density around 
the weakly interacting limit (fcpa)^^ — 0. 

Our results (see Figure [Sj 1/C is shown to be able to 
plot the repulsive and attractive side in the same picture) 
show good agreement with the exact relation ^. Eq. 
@ works well in the weakly to moderately interacting 
regime, although the agreement deteriorates at stronger 
interactions. Using the numerical derivative and eq. ([5]) 
gives better results. At k-pa = — 1 we obtain for the exact 
value of the contact Cox = (4 — 7r)/7r^ = 0.0870..., while 
using numerical values for the polaron energy and eq. 
dSl) gives Cnsc « 0.100 and eq. @ gives Cpert ^ 0.117. 
At k-pa = 1 we obtain Cox = (4 + 37r)/7r2 = 1.36..., 
Cnsc ~ 1-28 and Cpert ~ 1-57. No phase transition is 
found at any kpa] the jump in C at fcpa = is not a phase 
transition since the infinitely attractive and infinitely re- 
pulsive regime are not adiabatically connected. This is 
a peculiarity of the ID system; in the 3D case [30| there 
is a first order transition from a polaronic to a molecular 
state. The perturbative theory (j4]) assumes the polaron 
to be in the zero momentum state. Clearly, at stronger 
interactions, the polaron has an increasingly high prob- 
ability to be at higher momenta, and eq. (JH) becomes 
increasingly inaccurate. Correcting for this improves the 
results (squares in Figure [3]) and makes them approxi- 
mately in line with numerical results based on eq. ([5]) 
(stars in Figure [3]); details of this method are shown in 
the supplementary material. At fcpa = —1 this modified 
theory gives Cmod ~ 0.0970, while at kpa = 1 we obtain 



Cmod ~ 1-29. The difference between Cmod (squares) and 
Cnsc (stars) is due to an additional error from discretiz- 
ing the derivative in eq. ([5]) . The remaining error between 
Cmod a-nd Ccx is due to the NSC T-matrix approximation. 
At finite temperature (see Figure S]), the contact 
reaches the mean-field regime at high temperatures. In- 
terestingly, the contact reaches the same value regardless 
of the sign of kpCL in this limit. Therefore, physically, 
it does not matter whether one considers attractive or 
repulsive interactions in the high-temperature limit, as 
thermal excitations dominate. This feature can be un- 
derstood upon inspection of eq. (JH). As T — >■ oo, the 
density distribution will become spread out over higher 
momenta, and the value of the T-matrix at higher mo- 
menta, which corresponds to the mean-field behavior of 
the system, will become increasingly important. Another 
interesting feature is that for kpa < 0, the contact den- 
sity increases as the temperature increases. This is a ID 
feature that does not appear in the 3D spin-balanced case 
[9|. Once again, we find that as interactions increase, it 
takes higher temperatures to reach the mean-field regime. 
Physically, this means that while the mean-field predic- 
tion is a good approximation at high temperature for all 
values of fcpa, it takes increasingly high temperatures to 
reach this regime as \k-pa\ is reduced. 
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FIG. 4. (Color online) Contact density C as a function of in- 
verse temperature T-pjT . Top (green) line shows the attrac- 
tive polaron at fcpa = 2. Bottom (red) line shows the repulsive 
polaron at fcpa = —2. The limiting value at high temperature, 
where both lines meet asymptotically, is the same as the value 
for the contact indicated by the dotted line in Figure [S] Inset: 
normalized contact density C = r -r' for various values of 
kpa = —0.25, —0.5, —2, 2, —10, where "mf" is the mean-field 
prediction and "exact" is the contact density at T = for the 
considered scattering length. The bottom (cyan) line is for 
kpa = —0.25. All values were calculated using perturbation 
theory Q. 



In conclusion, we have computed numerical approxi- 
mations to the ID Fermi polaron ground state energies 



and contact densities over a large range of interaction 
strengths. In addition, we have derived an exact expres- 
sion for the contact density for the Fermi polaron in ID at 
zero temperature. Our numerical approximations show 
good agreement with exact results. Finally, we have com- 
puted numerical results for finite temperatures and show 
a crossover to the mean-field regime at high tempera- 
tures. 
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SUPPLEMENTARY MATERIAL TO THE LETTER: "Energy, contact and 
universal thermodynamics of the one-dimensional Fermi polaron at finite tempera- 
ture" 
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1 Contact from perturbation theory 

In Ref. [Ij one of the authors gave a perturbative result for calculating the prob- 
ability of finding the polaron in momentum state p (i.e. the polaron momentum 
distribution Vp) using the scattering T-matrix F. The one-dimensional analog of 
this result is 
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E 



Jlg{l -Uq^p), 



(1) 



where E is the polaron energy and rig is the majority component occupation number 
for momentum state q (we will here assume zero-temperature distributions and set 
?i = 1 for the sake of simplicity). The contact can be obtained from the high 
momentum limit of the momentum distribution as 



Cport = Im^ P'^'^p == TT y^ 



p— >oo 



E 



'2 



'(J) 



where in the non-self-consistent T-matrix theory 



r(go,'?) = 



9 



l-5E|fc' 



(2) 



(3) 



\k'\>kp 2qa-e^, 



Here g = —2/ma (this is the prefactor in front of the delta function in the inter- 
particle potential) and the sign convention implies that a is negative for repulsive 
interactions. As we shall soon see, this form for the contact is, up to a normalization, 
identical to the contact obtained from the first-order variational Chevy ansatz and 
in agreement with the Tan's adiabatic theorem. 



2 Contact from Chevy Ansatz 

The first-order variational Chevy Ansatz [5], in which only a single particle- hole 
excitation is considered, is 



|V)-0ofoS|O)+V0fe,6j_fccl.c,|O), 



kq 



(4) 



where the vacuum jO) corresponds to the ideal Fermi sea of majority atoms. The 
first term describes a zero-momentum polaron and the second term has the polaron 
(minority atom) in the momentum state q — k. The coefficients (f>kq will give the 
polaron momentum distribution as 



i^P = Wblbpli/j) = ^|0g_p,,|^ng(l -riq^p). 



(5) 



Solving the polaron momentum distribution is quite straightforward from the Chevy 
ansatz if one neglects the hole-hole scattering term, leading to equations f5] 

|fc|>fcF,|i3|<'i:F 

and 

94>0+9 X! (f'k'q = {E + eq - Ck- eq-k)4'kq, (7) 

|fc'|>/CF 

where n-f is the density of the majority component. From Eq. ([7]) we obtain 

**. . ^^;;;g';"'>f-' . (8) 

-Ci + Eg — e/c — Eg-fc 

Denoting the nominator in the right as x(g) — g4>Q + g '}2\k'\>k-f '^k'q yields 

which can be solved as 

500 J, -n f ^ + '^q (1 



xiq) = ^. ■" ^^ 1 = '^or ^-^, ^ . (10) 

The polaron momentum distribution can now be expressed by combining Eqs. ([S]), ([HI), 
and ((ini): 
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2 

n, (l-n^.p). (11) 



The expression on the right is the perturbative result in Eq. ^ except for the factor 
|(/)op. Indeed, the perturbative result is built on the assumption that the polaron is 
(with unit probability) in the zero-momentum state but the perturbative correction 
will yield an unnormalized state. Renormalizing the perturbed state would yield 
the correct prefactor and the exactly same result as in Eq. pT|) . Notice that by 
substituting Eq. pUj) into Eq. ^ would yield the polaron energy 

which is precisely the polaron energy obtained from the non-self-consistent T-matrix 
theory. 



3 Contact from energy 

The polaron energy, as computed from the Chevy ansatz, is 

E(a) = V ^^ ; = V ^ , . . (13) 



where g — —2/ via and energy E{a) will depend on scattering length a. According 
to the one-dimensional analog of the Tan adiabatic theorem as derived by Barth 
and Zwerger |4j 
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Calculating the derivative yields 
= ~-E{a) - ^ 
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[l-gx{E{a),q)]. 
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(15) 



where we have used the property gf = ^ f • Now 

d g dy 

^ [1 - 9X{E{a),q)] = -x{E{a),q) - 5-—, 
oa a oa 



which, upon substitution, yields 
dE 1 , , v^ 



9 



da 



q<kJ^-9XiEia),q)] 
Combining the first and second terms on the right yields 
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(18) 



The first term on the right is proportional to the perturbative result Cpcrt that we 
have been using in the main paper, the second term will yield the normalization as 
will be seen in a moment. Evaluating the differential on the second term yields 
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Thus we now obtain 
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q<fep [1 - 9XiEia), q)] ,^^,^^ [E{a) + e, - e^ - e,_fe 



r, (20) 



rearranging the sums in the second term and substituting the value of the contact 
yields 



c - Cpoit ~ ^ E E 

k'>kY q<.kF 

and further 

C = Cport - C' E E 
P#o q L 



9 



1 



[1 - 9X{E{a),q)]'^ [E{a) + e, - e^ - e<j-fe]^ 



71,(1 - nq^p)g'^ 



[1 - 9X{E{a), q)f [E{a) + e, - Cg-p - Cp]^ 



(21) 



(22) 



where the term in the square brackets is the perturbative result for the probabihty 
for a majority atom scattering from momentum state q to q — p. Oi course, should 
that be the case, the polaron will then have scattered to momentum state p and 
the summation over q will then give the polaron occupation number Vp. Thus we 
obtain 

C — Cpert — C 2_^ fpj (23) 

which finally yields 

C=—-^ . (24) 

That is, the perturbative result for the contact Cpcrt needs to be simply normalized 
by the polaron momentum distribution. In the weakly interacting limit, the proba- 
bility for the polaron not being in the zero momentum state is very small, and the 
perturbative result without the normalization correction is accurate. 
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